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Abstract 



C^ ' We study the response of a spin glass system with respect to the rescaling of 

its interaction random variables and investigate numerically the behaviour of the 
PJ ' correlation functions with respect to the volume. While for a ferromagnet the 

o: 

O . local energy correlation functions increase monotonically with the scale and, by 

consequence, with respect to the volume of the system we find that in a general 
spin glass model those monotonicities are violated. 

\o 
en 

^', 

O I When spins interact in a ferromagnetic system they tend to be aligned. That simple 

OO . 

"^ ■ fact is reflected in structural properties of the statistical mechanics equilibrium state called 



Griffiths inequalities |Gr[ IGr2] or more generally GKS |KS] inequalities. 

The first inequality states that the pressure of a spin system (minus the free energy 
times beta) does increase with the strength of each interaction among spins. The second 
says that the correlation among any set of spins increases with respect to the strength of 
the interaction of any other set of them. 

Since the strength of the interaction allows to switch on and off new parts of the 
system, the mentioned monotonicity properties can be easily turned into new ones like 
monotonicity with respect to the volume or with respect to the system dimensionality. 
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All those properties are at the origins of the fruitful applications of those inequalities to 
prove rigorous results in statistical mechanics, like the existence of the thermodynamic 
limit for pressure and correlations, bounds for critical temperatures and exponents [BSj 
and their mutual relations for different systems. 

While the physical meaning of the Griffiths inequalities was clear and well understood 
much earlier than rigorously proved, for disordered systems like spin glasses there are no a 
priori evident monotonicity properties due to the lack of ferromagnetism and the presence 
of competition (frustrated loops [MPVj ). 

The first case to be understood has been the Gaussian interaction. The use of the 
partial integration led to the first general proof of existence and monotonicity of thermo- 
dynamic limit for the pressure of a d-dimensional spin glass model with general potential 
|CG] . Defining the potential 

t/A(J,a) = 5^ Jxax, (0.1) 

where the coefficients Jx are a Gaussian family distributed as 

Av(Jx) = , Av(Jx^y) = A^5x,y (0.2) 

with its associated random Gibbs-Boltzmann state oj and the quenched measure as 

< - > = Av[cj(-)] (0.3) 

a straightforward computation jGGj gives 

< Jx(yx > = A^Av(l - uj,) , (0.4) 

which turns out to be positive by inspection. Although it was clear that the former 
positivity comes from convexity, it was recognized only few years later that the same 
result holds in full generality for random interaction Jx with zero average and not only 
for centered Gaussian variables: in ^CLj it is proved that (10. 4p can be simply derived from 
thermodynamic convexity. 

A natural perspective to look at the former inequalities is to consider the deformation 
of the general centered random variable Jx as XxJx with Xx > 0. The quenched pressure 
P as a function of the set of lambda's has first and second derivatives 



dXx 
d'^P d < Jx(Tx > 



dP 

< Jxcrx > > (0.5) 



r,. r,. o, Av[JxJy{^XY -LdxCdy)] (0.6) 

The two quantities (10. 5p and (10. 6p have been extensively studied in rf = 1 with nearest 
neighboor interaction and periodic (or free) boundary conditions in |CU] : the sign of (10.51) 
remains positive also by shifting on positive values the J averages; the value of (10.61) turns 
out to be non positive for zero mean and positive variance interactions and changes its 
sign when crossing a line in the mean-variance plane toward the ferromagnetic regime of 
zero variance and positive mean. 

The present paper deals with the study of the sign of (10. 6p for higher dimensions 
or different topologies for the case with zero average interaction. Our findings can be 
summarized as follows: (10. 60 doesn't have a definite sign. An explicit counterexample is 
found, for instance, in the case of a nearest neighboor spin chain with an extra interac- 
tion connecting two distant spins. In principle, a specific topology could not affect the 
monotonicity in the volume. For that reason we test numerically the nearest neighboor 
correlation function for two and three dimensional systems of increasing size and find an 
oscillating behaviour. 

Let consider a closed chain of six spins with nearest neighboor interaction with one 
added interaction between the spins 2 and 5. The dependence of the partition function 
on the couplings Ji^2 e J2,z is: 

a 

= acosh(/9Ji^2) cosh(/9J2,3) + 6sinh(/3Ji,2) cosh(/?J2,3) + 

+ccosh(/?Ji^2) sinh(/5J2,3) + (isinh(/?Ji_2) sinh(/?J2,3) 

where the four coefficients are: 



II cosh(/3J,J ■ 5^ n tanh(/3J,, 

Ji,j)^{l,2),{2,3) / 95=0 \{i,j)GB 



^{i,j)^{l,2),{2,3) J dB={l,2) \{i,j)€B 

c = [ n cosh(/5J,,,) I ■ 5^ I n tanh(/5J,,,) 

\ii,j)^il,2),{2,3) ) dB={2,3) \ii,j)€B 

d = i H cosh(/3J,,,) I ■ Yl ( n tanh(/5J,,,) 

\(i,j)7^{l>2),{2,3) / c»B=(l,2)U{2,3) \{i,j)€B 

The truncated correlation function is now: 

ac? — be 

(^12,23 — ^12 ^23 = 16 



Z2 

By Gauge invariance the Bernoulh random model can be reduced to one in which the 
randomness is concentrated on the two couplings Ji^2 = ilj <^2,3 = =tl with probability 
1/2, and the remaining others Jij = 1. 

Let introduce the notation: C := cosh(/9), S := sinh(/9), T := tanh(/3), and: T12 : = 
tanh(/9Ji,2), T23 := tanh(/9J2,3) (and similarly for sinh and cosh). 

We will indicate with Z{+, +), Z{+, — ), Z{—, +) and 2'(— , — ) the partition functions 
computed with fixed values of Ji 2 e J2,3. 
An explicit computation gives: 

and 

Z = C^ Ci C2 [1 + T12 T^ + T23 T3 + T12 T23 T^] => 

Z{++) = C^[l + 2T* + T^] 
Z{+-) = Z{-+) = C'il-T^] 

Z( ) = C^\1-2T* + T^] 



Finally: 



AV [71,2^2,3(^12,23 - ^12 ^0^23)] = 16 C' T^ (1 - T^) Av 



J\,2 J'. 



2,3 



Z2 



IQC'^'Tm-T' 



16^4 M _ rp2 



Z\+- 



C^ 



Z2(- 



+ 



Z2( 



1 + 2T4 + T6)2 (1-2T4 + T6)2 (1-T6)2 



16T4 (1 - T2) 8 T^ [-1 + 3T2 - 2T^ + 2TS - AT^^ - T 



3T 



141 



C4 (1 + 2T4 + T6)2 . (1 - 2T4 + T6)2 . (1 - T6)2 

In order to compute the sign, we notice that the square parenthesis term is: 



[-1 + 3T2 - 2T^ + 2T^ - 4Ti° - T^^ + 3T 



141 



^2(^2 - 1){C^ - 58) + 2S\C^ - S^) + S^iC^ + ^2) 
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c 

The only term with possible sign change is (52 — 1). For small (3 the leading term is then 
— cosh^°(/3), which gives a negative contribution, while for large P everything is positive. 
A plot of the function Av [712^2,3(1^12,23 — 1^12 (^23)] shows a change of sign around 
P = 0.695. 

The numerical test is performed for 7- dimensional cubic lattices A of volume A^ = L'^, 
with d = 2,3. We analyze two cases of quenched disorder: the Bernoulli couplings with 
Jj,j = ±1 and the Gaussian couplings with zero mean and unit variance. 

Given a spin configuration a for a system of linear size L, we consider the observable: 



Av[JbUJb] 



(0.7) 



where b = {i,j) with i,j & A, \i — j\ = 1, Ub is the thermal average of the quantity ajCTj 
and Av[-] is the average over the quenched disorder. 

With a parallel-tempering algorithm [HNj we investigate the correlation (10.71) for lat- 
tice sizes ranging from L = 3 to L = 24 in the case d = 2 and from L = 3toL = 10 
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Three dimensional lattice L"^, d = 3 
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Table 1: Parameters of the simulations: linear system size, number of sweeps used for 
thermalization, number of sweeps for measurement of the observable, number of disorder 
realizations, number of temperature values allowed in the parallel tempering procedure, 
temperature increment, minimum and maximum temperature values. 



in the case rf = 3. For each size we consider at least 2048 disorder realizations and, 
in order to thermalize the large sizes, we choose up to 37 temperature values in the 
range 0.5 < t < 2.3, in which the critical temperatures of the three dimensional models 
(tc — 0.95 for the Gaussian model |MPRLj and tc — 1.15 for the Bernoulh model |BCFj ) 
are contained. The thermalization in the parallel tempering procedure is tested by check- 
ing the symmetry of the probability distribution for the standard overlap q under the 
transformation q —>■ —q. Moreover, for the Gaussian coupling case it is available another 
thermalization test: the internal energy can be calculated both as the temporal mean of 
the Hamiltonian or, using integration by parts, as the expectation of a simple function of 
the link overlap [Coj . We checked that with our thermalization steps both measurements 
converge to the same value. All the parameters used in the simulations are reported in 
TabiH 

The numerical results are displayed in Fig {1] and in Fig 121 where the correlation (10.71) is 
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Figure 1: Af [JfcCjf,] as a function of L for the two dimensional lattice with gaussian and 
Bernoulli coupling, for two different temperature values t = 0.5 and t = 1.5 each. 

represented as a function of the linear system size L for different temperatures both for the 
two dimensional system and for the three dimensional one. We find that the correlation 
oscillates with respect to L, independently of the couplings (Bernoulli and gaussian), of 
the temperatures and of the dimension d {d = 2,3). 

The results presented here show that the monotonocity properties typical of ferro- 
magnetic systems are clearly violated for spin glass models. Further effort is necessary 
to establish if the quantity (10. 6p can keep a definite sign for specific lattice geometries or 
may depend on the relative position of the two sets X and Y. 
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